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DESCRIPTION 

TRANSFORMING A POLYGON I SED OBJECT INTO 
AN IMPLICIT -FUNCTION REPRESENTATION 

CROSS-REFERENCE TO RELATED APPLICATIONS 
This application is based upon and claims the 
benefit of priority from prior Japanese Patent 
Application No. 2004-231193, filed August 6, 2004, the 
entire contents of which are incorporated herein by 
reference . 

Technical Field 
The present invention relates generally to an 
image processing method, an image processing program 
product and an image processing apparatus, and more 
particularly to a scheme of transforming, for example, 
a polygonized object to an implicit-function represen- 
tation . 

Background Art 
In the prior art, in the field of computer 
graphics and CADs, an implicit-function representation 

w 

or a parametric representation using polygons, etc. has 
been used as means for representing the shape of an 
object. 

A method of modeling an organic object shape by 
an implicit-function representation is disclosed, 
for instance, in T. Nishita and E. Nakamae, "A Method 
for Displaying Metaballs by Using Bezie Clipping, 11 
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Proceeding of EUROGRAPHICS, 1994. In this method, a 
plurality of implicit surfaces, such as meatballs, 
are combined. In addition, a scheme of generating 
implicit surfaces by fitting curved surfaces to roughly 
5 defined point-group data is proposed in G. Turk and 

J.F. O'Brien, "Shape Transformation Using Variational 
Implicit Surface," Proceeding of SIGGRAPH, 1999. 

In these schemes using implicit surfaces, however, 
it is difficult to represent details of a complex 
10 character. If surface data that represents a detailed 

shape is to be generated, a high computational cost is 
required. 

In recent years, schemes called "level set 
methods", wherein implicit surfaces are generated 

15 using signed distance functions, have been proposed, 

for instance, in S. Osher and R. Fedkiw, "Level Set 
Methods: An Overview and Some Recent Results," Journal 
of Computational Physics, Vol. 169, 2001, and D. 
Enright, S. Marschner and R. Fedkiw, "Animation and 

20 Rendering of Complex Water Surface," Proceeding of 

SIGGRAPH, 2002. The advantage of using a signed 
distance function is that shape information can be 
represented with accuracy. Thus, when fluid animation 
is to be produced, this scheme is used as a represen- 

25 tation method for accurately handling a variation in 

surface shape of a fluid. 

The above prior art, however, merely discloses how 



to more accurately represent a variation in shape of 
an object that has a specified simple shape represented 
by an implicit function. The prior art is silent on a 
concrete scheme as to how to represent an object by an 
implicit function. There has been no proposal for a 
scheme for generally applying an implicit function to 
not only a simple-shape object but to a complex-shape 
ob j ect . 

Disclosure of Invention 
The present invention has been made in considera- 
tion of the above-described problems, and the object of 
the present invention is to provide an image processing 
method, a computer program product for processing image 
data and a image processing apparatus which can apply 
an implicit function to objects represented using 
polygons . 

An image processing method using a computer, 
according to an aspect of the present inven- 
tion, includes: extracting vertex coordinates of a 
triangular-shaped polygon; setting a region surrounding 
the triangular-shaped polygon on the basis of the 
vertex coordinates; measuring a distance from a lattice 
point included in the region to the triangular-shaped 
polygon; and drawing a graphic figure on the basis of 
the distance from the lattice point to the triangular- 
shaped polygon. 

A computer program product for processing image 



data, according to another aspect of the present 
invention, includes: means for instructing a computer 
to extract vertex coordinates of a triangular-shaped 
polygon; means for instructing the computer to generate 
a region surrounding the triangular-shaped polygon 
on the basis of the vertex coordinates; means for 
instructing the computer to measure a distance from a 
lattice point included in the region to the triangular- 
shaped polygon; and means for instructing the computer 
to draw a graphic figure on the basis of the distance 
from the lattice point to the triangular-shaped 
polygon . 

An image processing apparatus according to another 
aspect of the present invention includes: an input 
unit configured to receive polygon data; a processing 
unit configured to generate a region surrounding an 
individual polygon, which is represented by the polygon 
data, measure a distance from a lattice point included 
in the region to the polygon, and draw a graphic 
figure on the basis of the measured distance by an 
implicit-function representation; and an outputting 
unit configured to display the graphic figure that 
is obtained by the implicit-function representation. 

Brief Description of Drawings 

The file of this patent contains photographs 
executed in color. Copies of this patent with color 
photographs will be provided by the Patent and 
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Trademark Office upon request and payment of the 
necessary fee. 

FIG. 1 is a block diagram of an image processing 
computer according to a first embodiment of the present 
5 invention; 

FIG. 2 illustrates an image processing method 
according to the first embodiment of the invention, and 
is a conceptual view showing an object represented by 
polygons ("polygonal object") and an object represented 
10 by an implicit function ("implicit object"); 

FIG. 3 is a flow chart relating to the image 
processing method and an image processing program 
according to the first embodiment of the invention; 

FIG. 4 is a conceptual view illustrating a step 
15 of decomposing a rectangular-shaped polygon into 

triangular-shaped polygons in the image processing 
method and the image processing program according 
to the first embodiment; 

FIG. 5 is a conceptual view illustrating a step 
20 of obtaining vertex coordinates of a triangular-shaped 

polygon in the image processing method and the image 
processing program according to the first embodiment; 

FIG. 6 is a three-dimensional conceptual view 
illustrating a step of generating a bounding box in 
25 the image processing method and the image processing 

program according to the first embodiment; 

FIG. 7 is a two-dimensional conceptual view 



illustrating a step of generating the bounding box in 
the image processing method and the image processing 
program according to the first embodiment; 

FIG. 8 is a two-dimensional conceptual view 
illustrating a step of generating the bounding box in 
the image processing method and the image processing 
program according to the first embodiment; 

FIG. 9 is a two-dimensional conceptual view 
illustrating a step of generating a bounding sphere in 
the image processing method and the image processing 
program according to the first embodiment; 

FIG. 10 is a three-dimensional conceptual 
view illustrating a step of finding a signed minimum 
distance in the image processing method and the image 
processing program according to the first embodiment; 

FIG. 11 is a two-dimensional conceptual view 
illustrating a step of finding the signed minimum 
distance in the image processing method and the image 
processing program according to the first embodiment; 

FIG. 12 is a two-dimensional conceptual view 
illustrating a step of finding the signed minimum 
distance in the image processing method and the image 
processing program according to the first embodiment; 

FIG. 13 is a two-dimensional conceptual view 
illustrating a step of finding the signed minimum 
distance in the image processing method and the image 
processing program according to the first embodiment; 



FIG. 14 is a three-dimensional conceptual view 
illustrating a step of rendering an object by an 
implicit function in the image processing method and 
the image processing program according to the first 
embodiment ; 

FIG. 15 is a two-dimensional conceptual view 
illustrating a step of rendering the object by the 
implicit function in the image processing method and 
the image processing program according to the first 
embodiment ; 

FIG. 16 is a two-dimensional conceptual view 
illustrating a step of rendering the object by the 
implicit function in the image processing method and 
the image processing program according to the first 
embodiment ; 

FIG. 17 is an image photograph of an object on 
a display, which is rendered by the image processing 
method and the image processing program according to 
the first embodiment; 

FIG. 18 is an image photograph of an object on 
the display, which is rendered by the image processing 
method and the image processing program according to 
the first embodiment; 

FIG. 19 is an image photograph of an object on 
a display, which is rendered by a prior-art polygonal 
representation; 

FIG. 20 is an image photograph of an object on the 
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display, which is rendered by the prior-art polygonal 
representation; 

FIG. 21 is an image photograph of a cross section 
in an X-Y plane of an object on a display, which is 
5 rendered by the image processing method and the image 

processing program according to the first embodiment; 

FIG. 22 is an image photograph of a cross section 
in an X-Z plane of the object on the display, which is 
rendered by the image processing method and the image 

10 processing program according to the first embodiment; 

FIG. 23 is an image photograph of a cross section 
in a Y-Z plane of the object on the display, which is 
rendered by the image processing method and the image 
processing program according to the first embodiment; 

15 FIG. 24 is a two-dimensional conceptual view 

illustrating a state in which a bounding box is 
generated in an image processing method and an image 
processing program according to a second embodiment 
of the invention; 

20 FIG. 25 is a flow chart relating to the image 

processing method and the image processing program 
according to the second embodiment of the invention; 

FIG. 26 is a two-dimensional conceptual view 
showing an object that is transformed from a polygonal 

25 representation to an implicit-function representation 

by the image processing method and the image processing 
program according to the second embodiment; 



« 
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FIG. 27 is a two-dimensional conceptual view 
illustrating a step of generating a bounding box in 
the image processing method and the image processing 
program according to the second embodiment; 
5 FIG. 28 is a two-dimensional conceptual view 

illustrating a step of dividing a bounding box in 
the image processing method and the image processing 
program according to the second embodiment; 

FIG. 29 is a two-dimensional conceptual view 
10 illustrating a step of dividing a bounding box in 

the image processing method and the image processing 
program according to the second embodiment; 

FIG. 30 is a two-dimensional conceptual view 
illustrating a step of dividing a bounding box in 
15 the image processing method and the image processing 

program according to the second embodiment; 

FIG. 31 is a two-dimensional conceptual view 
illustrating a step of dividing a bounding box in 
the image processing method and the image processing 
20 program according to the second embodiment; 

FIG. 32 is a flow chart relating to a image proc- 
essing method and a image processing program according 
to a third embodiment of the invention; 

FIG. 33 is a three-dimensional conceptual 
25 view showing a bounding box that is generated by the 

image processing method and image processing program 
according to the third embodiment, illustrating a state 
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in which a plurality of polygons are included in one 
bounding box; 

FIG. 34 is a two-dimensional conceptual view 
illustrating a step of finding a signed minimum 
distance in the image processing method and the image 
processing program according to the third embodiment, 
FIG. 34 being an X-Y plane at a time when the signed 
minimum distance is found using two polygons; 

FIG. 35 is a two-dimensional conceptual view 
illustrating a step of finding a signed minimum 
distance in the image processing method and the image 
processing program according to the third embodiment, 
FIG. 35 being an X-Y plane at a time when the signed 
minimum distance is found using two polygons; 

FIG. 36 is a block diagram of an image proc- 
essor system LSI that includes an image processing 
method, an image processing program and an image 
processing computer according to a fourth embodiment 
of the invention; 

FIG. 37 is a block diagram of an arithmetic 
process unit in the image processor system LSI that 
includes the image processing method, image processing 
program and image processing computer according to the 
fourth embodiment of the invention; 

FIG. 38 is a block diagram of a digital board of a 
digital TV that includes the image processing method, 
image processing program and image processing computer 



according to the first to third embodiments; and 

FIG. 39 is a block diagram of a recording/ 
reproducing apparatus that includes the image 
processing method, image processing program and image 
processing computer according to the first to third 
embodiments . 

Best Mode for Carrying Out the Invention 
Embodiments of the present invention will now be 
described with reference to the accompanying drawings. 
In the embodiments described below, an arrow used in, 
e.g. " (r— *)", represents a vector, and does not mean 
"derivation" in a symbol logic or "image" in a set 
theory . 

Referring to FIG. 1, a description is given of an 
image processing method, an image processing program 
product and an image processing computer according to a 
first embodiment of the invention. FIG. 1 is a block 
diagram showing the internal configuration of the image 
processing computer according to the embodiment. 

As is shown in FIG. 1, an image processing 
computer 10 according to the embodiment comprises a 
CPU 11, a main memory 12, a memory 13, an input device 
14, an output device 15 and a data bus 16. 

When an image process is executed, the CPU 11 
reads out an image process program 17 from the memory 
13 and loads it in the main memory 12. The CPU 11 
executes the image process according to the image 



process program 17 that is loaded in the main memory 
12. In the image process, the main memory 12 stores 
data, which is frequently used by the CPU 11, as well 
as the image process program product 17. The memory 13 
is, e.g. a readable/writable nonvolatile semiconductor 
memory. The memory 13 stores various programs and 
data. The input device 14 receives data from outside. 
The output device 15 is, e.g. a display or a printer, 
which outputs an image that is rendered by the CPU 11. 
The data bus 16 executes transmission/reception of 
data between respective blocks in the image processing 
computer 10. 

Next, the operation of the image processing 
computer 10 with the above configuration is described. 
The image processing computer 10 transforms a para- 
metrically represented object, such as a polygonized 
object, to an implicit-function representation. FIG. 2 
illustrates the scheme of this conversion. FIG. 2 
schematically shows an object represented by polygons 
("polygonal object") and an object represented by an 
implicit function ("implicit object") . Assume now that 
a spherical object is represented by triangular-shaped 
polygons, as shown in FIG. 2. The image processing 
computer 10 represents the spherical object using an 
implicit function . 

The implicit-function representation is the 
following scheme. In a case where a surface S C r3 of 



an object is given, the implicit-function representa- 
tion of the surface S is a set {t, d>(r— >)}, which meets 

r-* e S <=> <J> (r->) = t 
where t is a given real value, and <t> (r— >) is a real- 
valued function. In this case, r— * e R 3 represents a 
position in a three-dimensional space. Without losing 
generality, t = 0 may be set. Normally, the surface 
of an object is defined by 0 (r— *) =0. In the present 
embodiment, <f> (r— >) is used not in the form of a 
function, but in the form of a surface position of an 
object. In this case, the surface position of the 
object is expressed using a signed minimum distance. 
The function using the signed minimum distance provides 
a value with a plus sign when a nearest point is 
present on the obverse side of a surface at a distance 
from a given point r— > to the nearest point on the 
surface S, and provides a value with a minus sign when 
the nearest point is present on the reverse side of the 
surface. The signed minimum distance will be described 
later in greater detail. 

A concrete method of transforming a polygonal 
representation to an implicit-function represen- 
tation using the image processing computer 10 is 
described referring to FIG. 3. FIG. 3 is a flow chart 
illustrating a method of transforming a polygonal 
representation to an implicit-function representation 
using the image processing computer according to the 



present embodiment . 

As is shown in FIG. 3, to start with, data on a 
polygonized object, or an polygonal object, is input 
through the input means 14 shown in FIG. 1 and stored, 
for example, in the memory 13 of the image processing 
computer 10 (step S10) . 

Then, the CPU 11 reads out the image process 
program 17 from the memory 13 and loads it in the main 
memory 12. The image process program 17 may be read 
out from a portable magnetic memory medium, instead of 
the memory 13. Further, the image process program 17 
may be downloaded from a Web site via the Internet. 
At the same time, the CPU 11 reads out the object data 

(that is stored in step S10) from the memory 13 and 
loads it in the main memory 12. The CPU 11 executes an 
image process for the polygonal object according to the 
image process program 17 that is loaded in the main 
memory 12, thereby to transform the polygonal object 
to an implicit-function representation. A series of 
process steps of the image process are described below. 

The CPU 11 confirms, based on the object data 
loaded in the main memory 12, whether all the polygons 
of the object surface are triangular-shape polygons 

(step Sll) . If there is a polygon with four or more 
angles (n-angular polygon; n > 3) , the CPU 11 divides 
the polygon into triangular-shaped polygons (step S12) . 
FIG. 4 illustrates this process step. FIG. 4 



illustrates a case where a rectangular-shaped polygon 
is present on the surface of a spherical object. 
The CPU 11 divides the rectangular-shaped polygon into 
two triangular-shaped polygons. Needless to say, the 
same process is executed in a case where a pentagonal- 
shaped polygon is included or in a case where a 
hexagonal-shaped polygon is included. In these cases, 
the polygons are divided into three triangular-shaped 
polygons and four triangular-shaped polygons, respec- 
tively. By the process in" step S12, all the polygons 
that form the object surface are made into 
triangular-shaped polygons. 

Next, the CPU 11 extracts vertex coordinates of 
each triangular-shaped polygon (step S13) . FIG. 5 
illustrates this process step. As shown in FIG. 5, 
assume that there is a triangular-shaped polygon TP 
having three vertices PI, P2 and P3 that are expressed 
by rl— >, r2— > and r3— ► from the origin 0. The CPU 
12 finds coordinates (xl, yl, zl), (x2, y2, z2) and 
(x3, y3, z3) of vertices PI, P2 and P3 in a three- 
dimensional space. This process is executed for all 
the triangular-shaped polygons that represent the 
object surface. 

Subsequently, based on the vertex coordinates 
found in step S13, the CPU 11 generates a bounding box 
BB or a bounding sphere BS (step S14) . FIG. 6 three- 
dimensionally shows an example of the bounding box BB . 



As is shown in FIG. 6, the bounding box BB refers to 
a space that surrounds one triangular-shaped polygon 
TP. When the bounding box BB or bounding sphere BS 
is generated, the CPU 11 generates lattice points that 
are arranged at regular intervals. In the case of 
the bounding box BB, it is conceptually assumed that 
lattice points LP define a space, i.e. a cube, that 
surrounds the triangular-shaped polygon. The lattice 
points LP are arranged at least at the corners of the 
cube. The bounding box BB or bounding sphere BS is 
generated as a set of the lattice points. FIG. 7 two- 
dimensionally shows an example of the bounding box BB. 
As shown in FIG. 7, a plurality of lattice points LP 
are arranged in a matrix in a three-dimensional space 
in which the triangular-shaped polygon TP is present. 
Since FIG. 7 depicts the space two-dimensionally , 
the triangular-shaped polygon can be viewed as a line 
segment . 

As is shown in FIG. 8, the bounding box BB is 
generated with a size that is greater than a size of 
the triangular-shaped polygon TP by a degree corre- 
sponding to some lattice points ( e ) . FIG. 9 shows a 
case of the bounding sphere BS. Like the bounding 
box BB, the bounding sphere BB is generated so as to 
surround the triangular-shaped polygon TP with a size 
that is greater than the size of the triangular-shaped 
polygon TP by the degree ( a ) . In the description 



below, the bounding box BB is described by way of 
example, but the same applies exactly to the case of 
the bounding sphere BS. 

Next, the CPU 11 measures a signed minimum 
distance from the lattice point LP included in the 
bounding box BB to the triangular-shaped polygon TP 
included in the bounding box BB (step S15) . The method 
of measuring the signed minimum distance D(r— >) is 
described referring to FIG. 10 to FIG . 13. FIG. 10 
is a schematic diagram that three-dimensionally shows 
a plane including the plane of the triangular-shaped 
polygon TP. The plane including the plane of the 
triangular-shaped polygon TP is defined as an X-Y 
plane, and a normal vector to the X-Y plane is defined 
as a Z axis. FIG. 11 and FIG. 13 show X-Z planes, and 
FIG. 12 shows an X-Y plane. In FIG. 11 and FIG. 13, a 
direction extending from an intersection between the X 
axis and Z axis perpendicularly to the sheet surface 
of the drawings is a Y axis direction. In FIG. 12, a 
direction extending from an intersection between the X 
axis and Y axis perpendicularly to the sheet surface of 
the drawing is the Z axis direction. 

As is shown in FIG. 10 and FIG. 11, when a lattice 
point Al or A2 is arbitrarily given, this lattice 
point is set as r— > = (x, y, z) . The origin of r— » is 
set as an origin O. A nearest point of r— », which is 
present on the triangular-shaped polygon TP, is set as 
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r n — » = (x n/ y n , z n ) . Then, as shown in FIG. 10, any 
one of the vertices of the triangular-shaped polygon TP 
is set as an origin O 1 , and r— > is transformed to r 1 —> = 

(x ! , y ! , z') on a coordinate system having an X axis, 
which extends from the origin O' and corresponds to one 
of the sides of the triangular-shaped polygon TP, a Y 
axis, which is perpendicular to the X axis in the plane 
including the triangular-shaped polygon TP, and a Z 
axis that corresponds to the normal vector to the plane 
including the triangular-shaped polygon TP. Further, a 
projection point of r 1 — > onto the X-Y plane is set as 
ro f ^ = (x 1 , y ! , 0). In the triangular-shaped polygon 
TP, a nearest lattice point from rg 1- * becomes r n — > . 
As is clear from FIG. 11, the absolute value of the 
distance from the arbitrary lattice point to the 
nearest lattice point is given by I r— ► — r n — > | . 
A signed minimum distance D(r— >) is given by sgn(z ! ) 

| r—* — r n — > | . In this case, sgn(z f ) is a function 
indicative of whether the value z 1 of the lattice point 
of interest is plus or minus on the XYZ coordinates. 

Taking the lattice points Al and A2 by way of 
example, a specific method of measuring the signed 
minimum distance D(r— *) is explained. The case of the 
lattice point Al is referred to as CASE 1, and the case 
of the lattice point A2 as CASE 2. CASE 1 is first 
described . 

To begin with, the position r— > of the lattice 
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point Al is transformed to r 1 — » on the XYZ coordinate 
system. The lattice point Al is projected onto the X-Y 
plane. The projection position is ro f ~~> = (x f , y f , 0) . 
In CASE 1, as shown in FIG. 10 and FIG. 11, r 0 '^ is 
5 present within the plane of the triangular-shaped 

polygon TP. Thus, r 0 '^ = (X', y', 0) = r n — ► = (x n , y n , 
z n ) , and the position r 1 — > of the lattice point Al 
coincides with the position r n — * of a nearest lattice 
point Bl on the X-coordinate and Y-coordinate . Thus, 
10 the signed minimum distance D(r— >) at the lattice 

point Al becomes sgn(z ! ) | r— ► — r n~ * I = + zf (see 
FIG. 11) . 

Next, CASE 2 is described. Like CASE 1, the 
position r— » of the lattice point A2 is transformed to 

15 r'— > on the XYZ coordinate system. The lattice point 
A2 is projected onto the X-Y plane. The projection 
position is rQ ,_> = (x 1 , y f , 0). In CASE 2, as shown 
in FIG. 10 and FIG. 11, rg 1 — ► is present outside 
the triangular-shaped polygon TP. Thus, the position 

20 r n — * = ( x n* Yn^ z n) °f a nearest point B2 associated 
with the lattice point A2 is different from the 
projection position rg' - > = (x 1 , y ! , 0) of the 
lattice point A2 . The nearest point B2 is located 
at a point to which the projection point rQ l— > is 

25 shifted perpendicular to one side of the triangular- 

shaped polygon TP (see FIG. 10) . Thus, the signed 
minimum distance D(r— at the lattice point A2 is 
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sgn(z') | r— > — r n — » I = + I r— > — r n ^ I z') (see 

FIG . 11) . 

As has been described above, in CASE 1, the 
nearest lattice point is present within the plane of 
the triangular-shaped polygon TP. On the other hand, 
in CASE 2, the nearest lattice point is present on 
the side of the triangular-shaped polygon TP. As 
regards this point, a description is given referring to 
FIG. 12. Taking into account the relationship between 
an arbitrary lattice point and a lattice point nearest 
to the arbitrary lattice point on the triangular-shaped 
polygon, the X-Y plane can be divided into three areas 
AREA1, AREA2 and AREA3, as shown in FIG. 12. The area 
AREA1 is an area within the plane of the triangular- 
shaped polygon TP. The area AREA2 is an area defined 
between two parallel straight lines that extend from 
the associated vertices of the triangular-shaped 
polygon TP perpendicularly to the associated sides 

■ 

thereof. The area AREA3 is an area other than the 
areas AREA1 and AREA2 . 

When an arbitrary lattice point is projected on 
the X-Y plane, if the projection position is present 
within the area AREA1, ro'"* = r n — *, and r n — > is 
present within the plane of the triangular-shaped 
polygon TP. This case corresponds to the above- 
described CASE 1. 

If the projection position is present in the 



area AREA2, ro'— > r n ^, and r n ~^ is present on 
the side of the triangular-shaped polygon TP. A 
displacement between rQ ,— > and r n — ► in the X-Y plane 
is ( | r n — » — rg 1 ^! 2 ) 1 / 2 . This case corresponds to 
the above-described CASE 2. 

Further, if the projection position is present in 

the area AREA3, rg 1- ^ r n~*' and r n^ is present at 
the vertex of the triangular-shaped polygon TP. A 
displacement between rg 1- * and r n — > in the X-Y plane 

is ( |r n -> - r 0 '^l 2 ) 1/2 . 

In the example shown in FIG. 11 and FIG. 12, the 
Z-axis value of the lattice point Al , A2 is plus. In 
FIG. 13, the Z-axis value of the lattice point Al is 
minus. In this case, sgn(z f ) is minus. The plus/minus 
of the signed function sgn(z') indicates whether the 
lattice point is present on the outside or on the 
inside of the triangular-shaped polygon TP. Specifi- 
cally, if the sgn(z') is plus, the lattice point Al 
is present on the outside of the object surface 
represented by the triangular-shaped polygon TP. If 
the sgn(z ? ) is minus, the lattice point Al is present 
on the inside of the object surface represented by the 
triangular-shaped polygon TP. 

The process of step S15 is executed, as described 
above. The process of step S15 is repeated for all the 
lattice points LP that are included in the bounding box 
BB (step S16, S17) . When the process of step S15 is 



completed for all lattice points LP, the signed minimum 
distance D(r— ►) associated with the bounding box BB is 
determined (step S18) . 

As has been described above, in steps S14 to S18, 
the signed minimum distance D(r— ►) associated with one 
triangular-shaped polygon TP is found. The CPU 11 
repeats execution of the process of steps S14 to 
S18 for all the triangular-shaped polygons (steps 
S19, S20) . If the process for all triangular-shaped 
polygons is completed, the CPU 11 renders the to-be- 
processed object on the basis of the signed minimum 
distance D(r— ») associated with each triangular-shaped 
polygon (step S21) . 

The object is rendered as a set of points with the 
signed minimum distances D(r— *) =0. In step S15, the 
distance between each lattice point and the triangular- 
shaped polygon has been found. Thus, by overlapping 
the distances from the respective lattice points 
to the triangular-shaped polygon, the position of 
the surface of the triangular-shaped polygon can be 
found. In short, the set of points with the signed 
minimum distances D(r— ») = 0 does meet the implicit 
function <f> (r— *) = 0, and does become the solution of 
the implicit function that represents the surface S of 
the to-be-processed object. Hence, in the present 
embodiment, the implicit function $ (r— ») itself 
is not found, but the set of signed minimum distances 



D(r— ►) = 0 is found, thereby substantially transforming 
the object to the implicit function <t> (r— >) . 

As described above, the signed minimum distances 
D(r— ») is found for each polygon. Consequently, there 
may be a case where the relation between adjacent 
polygons cannot be maintained by merely overlapping the 
signed minimum distances D(r— FIG. 14 and FIG. 15 
illustrate such a case. FIG. 14 and FIG. 15 are 
schematic views illustrating the state in which signed 
minimum distances D(r— ») are found for adjacent 
triangular-shaped polygons TP1 and TP2. FIG. 14 is a 
three-dimensional schematic view, and FIG. 15 is a 
two-dimensional schematic view. 

Assume, as shown in FIG. 14 and FIG. 15, that the 
two triangular-shaped polygons TP1 and TP2 are adjacent 
to each other at their vertices PI and P2 . However, as 
shown in FIG. 14, the signed minimum distances D(r— ►) 
associated with the triangular-shaped polygons TP1 and 
TP2 are found by different bounding boxes BB1 and BB2 . 
As a result, as shown in FIG. 15, there is a possi- 
bility that the vertices PI and P2 do not agree when 
the polygonal representation is transformed to the 
implicit-function representation. In this case, as 
shown in FIG. 16, the CPU 11 adjusts the positions of 
the object surface so that the vertices PI and P2 may 
agree. In FIG. 15 and FIG. 16, hatched areas are the 
areas where the signed function sgn(z f ) is minus. 



If the two polygons TPl and TP2 are made to agree, as 
shown in FIG. 16, the area where the signed function 
sgn(z') is minus is confined in the region surrounded 
by the two polygons TPl and TP2. In other words, when 
the spherical object, for example, as shown in FIG. 2 
is considered, the signed minimum distances D(r— ») = 0 
are set on the surface S of the spherical object, the 
area where the signed function sgn(z ! ) is minus is 
confined in the area on the inside of the spherical 
object, and the area where the signed function sgn(z ! ) 
is plus is present in the area on the outside of the 
spherical object . 

By the above process, the polygonal object is 
transformed to the implicit-function representation, 
and the CPU 11 outputs the object, which is represented 
by the implicit function, through the output device 
15. FIG. 17 and FIG. 18 show examples of objects 
that are represented by the implicit functions by the 
method described in the present embodiment. The 
transformation to the implicit function is performed 
using (90 X 90 X 90) bounding boxes. FIG. 19 and 
FIG. 20 show examples of objects that are obtained by 
polygonizing the objects shown in FIG. 17 and FIG. 18. 
Each of these polygonal objects is represented by 8754 
vertices and 17504 surfaces. As is clear from the 
comparison between FIGS. 17 and 18 and FIGS. 19 and 
20, the object surface is smoother and more natural in 



25 



the implicit function representation of the present 
embodiment than in the polygonal representation. 

FIG. 21 through FIG. 23 are an X-Y cross section, 
an X-Z cross section and a Y-Z cross section of the 
5 object shown in FIG. 17 and FIG. 18. In the Figures, 

an area colored in red is an area where the signed 
function sgn(z') is minus, an area colored in blue 
is an area where the signed function sgn(z') is plus, 
and an area colored in green is an area where the 

10 signed function sgn(z') is nearly zero. As shown in 

the Figures, the area with the minus-signed function 
sgn(z f ) is confined inside the object, and the area 
with the plus-signed function sgn(z') is present on the 
outside of the object. In this manner, a surface which 

15 has, in fact, no thickness may be made to have a 

thickness, and thus a phenomenon that occurs near the 
object surface may be treated in a simplified manner. 

The following advantageous effects (1) and (2) 
can be obtained by the image processing method, image 

20 processing program product and image processing 

computer according to the present embodiment. 

(1) An implicit-function representation method 
with high versatility in use can be obtained. 

In the method of this embodiment, the signed 

25 distance function D(r— ») is found for each of indi- 

vidual triangular-shaped polygons TP. Using a set 
of points with D(r— ») = 0, the object is represented by 
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the implicit function <t> (r— >) . Accordingly, regardless 
of the shape of the object to be rendered, the 
polygonal representation can be transformed to 
the implicit-function representation. In other words, 
5 no matter how complex the shape of the object is, 

the implicit-function representation method with high 
versatility in use can be provided. 

In addition, since the implicit-function repre- 
sentation can be used with high versatility, animation 
10 with more reality can be created in various fields, 

compared to the polygonal representation. 

(2) Determination of collision of objects can be 
simplified. 

The implicit-function representation according to 
15 this embodiment uses the signed distance function. The 

signed function sgn(z') is minus on the inside of the 
object to be rendered, and it is plus on the outside of 
the object. Thus, determination of collision between 
two objects is very simple. For example, if an object 
20 B collides with an object A, the surface of the object 

B contacts the object A, or a part of the surface of 
the object B is present inside the object A. This can 
easily be determined by paying attention to the signed 
function sgn(z') relating to the surface of the object 
25 B. Therefore, computer games, such as racing games or 

fighting games, that includes more complex objects, can 
be enjoyed. 



Next, a description is given of an image 
processing method, an image processing program product 
and an image processing computer according to a 
second embodiment of the invention. In this second 
embodiment, bounding boxes are set more finely in a 
region including an object outline part than in the 
first embodiment. 

FIG. 2 4 is a schematic diagram that two- 
dimensionally illustrates a state in which bounding 
boxes BB are generated with respect to a polygonal 
object OB. Squares in FIG. 24 are bounding boxes BB. 
As is shown in FIG. 24, in the present embodiment, 
the size of the bounding box is adjusted on a location- 
by-location basis. In particular, the size of the 
bounding box is decreased at an outline part of the 
object OB to be treated, and a greater number of 
bounding boxes BB are generated. 

Referring to a flow chart of FIG. 25, how to find 
a signed minimum distance D(r~ >) in this embodiment is 
described. The image processing computer according to 
the present embodiment has the same structure as that 
of the first embodiment, so a description thereof is 
omitted. To start with, the process up to step Sll or 
S12 described in connection with the first embodiment 
is executed, and all the polygons included in the to- 
be-treated object OB are divided into triangular-shaped 
polygons. Assume that the to-be-treated object OB 
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has a shape as shown in, e.g. FIG. 26. FIG. 26 
two-dimensionally shows the object OB. 

Subsequently, according to the image process 
program 17 that is loaded into the main memory 12, 
the CPU 11 resets a lattice-setting number-of -times k 
(step S30) . The lattice-setting number-of -times k is a 
numerical value indicating how finely the bounding box 
is to be divided, as will be described later in detail. 
The upper limit value of the lattice-setting number- 
of-times k is preset, and the CPU 11 causes the main 
memory 12 to store the lattice-setting number-of -times 
k and the upper limit value thereof. 

As is shown in FIG. 27, the CPU 11 divides a 
region, which includes the object OB, into a plurality 
of meshes Ml to M25. Each mesh is a cube having a 
length L on each side, A bounding box BB is generated 
for each mesh (step S31) . Lattice points included in 
the bounding box BB may be those as described in the 
first embodiment, or may be generated only at the eight 
vertices of the cube. FIG. 27, too, two-dimensionally 
shows the object OB. If the region that includes the 
to-be-treated object OB is a cube, cubes are arranged 
in five columns in a direction perpendicular to the 
sheet surface of FIG. 27. The total number of gener- 
ated bounding boxes BB is (25 X 5 columns) = 125. For 
the purpose of simple description, in the following 
descriptions, attention is paid to only 25 meshes Ml to 
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M25 that appear in the two-dimensional space . 

The CPU 11 successively finds the signed minimum 
distance D(r— >) from a lattice point to a triangular- 
shaped polygon with respect to the bounding boxes BB 
5 corresponding to the meshes Ml to M25 (steps S15 to 

S17). How to find the signed minimum distance D(r— *) 
has been described in connection with the first 
embodiment, so a description thereof is omitted here. 
If the signed minimum distances D(r— ») associated 

10 with all lattice points in an arbitrary bounding box BB 

(step S16) , the CPU 11 discriminates whether the signed 
functions sgn(z') associated with all lattice points are 
plus or not (step S32) . If all signed functions sgn(z') 
are plus, the signed minimum distances D(r— ») in this 

15 bounding box BB are determined. For the next bounding 

box BB (step S39) , the process of steps S15 to S17 is 
executed. If it is discriminated in step S32 that all 
signed functions sgn(z ? ) are not plus, the CPU 11 then 
discriminates whether all the signed functions sgn(z f ) 

20 are minus or not (step S33) . If all the signed 

functions sgn(z f ) are minus, the signed minimum 
distances D(r—>) in the bounding box BB are determined. 
For the next bounding box BB (step S39) , the process of 
steps S15 to S17 is executed. In the example shown in 

25 FIG. 27, if the process starts with the mesh Ml, the 

signed minimum distances D(r— >) for the meshes Ml to M6 
are determined by the above process. 
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If it is determined in step S33 that all the 
signed functions sgn(z') are not minus, that is, if 
the bounding box BB includes lattice points with 
plus-signed functions sgn(z f ) and lattice points with 
minus-signed functions sgn(z'), the bounding box BB 
includes an outline part of the object (step S34) . 
In the example of FIG. 27, this case first occurs with 
respect to the bounding box BB associated with the mesh 
M7 . Then, the CPU 11 compares the lattice-setting 
number-of -times k with the upper limit value thereof 
that is stored in the main memory 12 (step S18). If 
the upper limit value is reached, the latest found 
values are determined as the signed minimum distances 
D(r— *) (step S18) . If the upper limit value is not 
reached, the CPU 11 divides the mesh corresponding to 
this bounding box BB into finer meshes each having a 
1/2 side dimension (step S36) . FIG . 28 illustrates 
this process step. FIG. 28 illustrates the process 
of step S36 for the mesh M7 . The mesh M7 is divided 
into meshes M7-10 to M7-13 each having an L/2 side 
dimension. The CPU 11 generates bounding boxes BB in 
association with the meshes M7-10 to M7-13 (step S36) , 
and increments the lattice-setting number-of -times k by 
+1. The lattice-setting number-of -times k indicates 
the number of times of re-division of a mesh, and the 
upper limit value of the lattice-setting number-of- 
times k determines how finely the mesh is set. In 
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this description, it is assumed that the upper limit 
value is 3. 

The CPU 11 repeats the process of steps S15 to S35 
for the bounding boxes BB associated with the meshes 
5 M7-10 to M7-13, into which the mesh M7 is divided. As 

shown in FIG. 28, an outline part of the object OB is 
included in only the mesh M7-13 in the mesh M7 . 
Accordingly, all the signed functions sgn(z') of the 
lattice points in the bounding boxes BB corresponding 

10 to the meshes M7-10 to M7-12 are plus (step S32), and 

thus the signed minimum distances D(r— > ) for these 
bounding boxes BB are determined (step S18) . On the 
other hand, the bounding box BB corresponding to the 
mesh M7-13 includes both lattice points with plus- 

15 signed functions sgn(z ! ) and minus-signed functions 

sgn(z') (step S32, S33). Thus, the CPU 11 re-divides 
the mesh M7-13 into meshes M7-20 to M7-23 each having 
an L/4 side dimension, and generates bounding boxes 
BB corresponding to these meshes M7-20 to M7-23 (step 

20 S36) . FIG. 29 illustrates this process. The lattice- 

setting number-of -times k is incremented to k = 2 
(step S37 ) . 

Then, the CPU 11 repeats the process of steps S15 
to S35 for the bounding boxes BB associated with the 
25 meshes M7-20 to M7-23, into which the mesh M7-13 is 

re-divided. As shown in FIG. 29, an outline part of 
the object OB is included in only the mesh M7-23 in 
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the mesh M7-13. Accordingly, all the signed functions 
sgn(z') of the lattice points in the bounding boxes 
BB corresponding to the meshes M7-20 to M7-22 are 
plus (step S32), and thus the signed minimum distances 
D(r— ») for these bounding boxes BB are determined (step 
S18) . On the other hand, the bounding box BB 
corresponding to the mesh M7-23 includes both lattice 
points with plus-signed functions sgn(z') and minus- 
signed functions sgn(z ! ) (step S32," S33) . Thus, the 
CPU 11 re-divides the mesh M7-23 into meshes M7-30 to 
M7-33 each having an L/8 side dimension, and generates 
bounding boxes BB corresponding to these meshes M7-30 
to M7-33 (step S36) . FIG . 30 illustrates this process. 
The lattice-setting number-of -times k is incremented to 
k = 3 (step S37) . 

Subsequently, the CPU 11 repeats the process of 
steps S15 to S35 for the bounding boxes BB associated 
with the meshes M7-30 to M7-33, into which the mesh 
M7-23 is re-divided. As shown in FIG. 30, an outline 
part of the object OB is included in only the mesh 
M7-33 in the mesh M7-23. Accordingly, all the signed 
functions sgn(z f ) of the lattice points in the bounding 
boxes BB corresponding to the meshes M7-30 to M7-32 are 
plus (step S32), and thus the signed minimum distances 
D(r— >) for these bounding boxes BB are determined 
(step S18) . On the other hand, the bounding box BB 
corresponding to the mesh M7-33 includes both lattice 
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points with plus-signed functions sgn(z') and minus- 
signed functions sgn(z') (step S32, S33). In this 
case, however, the lattice-setting number-of -times k 
already reaches the upper limit value of k = 3 (step 
5 S35) . Thus, the CPU 11 also determines the signed 

minimum distances D(r— *) for the bounding box BB 
corresponding to the mesh M7-33 (step S18). 

The process for the mesh M7 is thus completed, and 
the CPU 11 executes the process of steps S15 to S37 for 

10 the bounding box BB corresponding to the mesh M8 (step 

S38, S39) . Since the mesh M8, too, includes an outline 
part of the object OB, the CPU 11 repeats, like the 
case of the mesh M7, re-division of associated divided 
meshes until the lattice-setting number-of -times k 

15 reaches the upper limit value. Thereby, the CPU 11 

determines the signed minimum distances D(r— 
FIG. 31 illustrates this process. 

Subsequently, similar processes are repeated for 
the meshes M9 to M25. As is shown in FIG. 27, all the 

20 signed functions sgn(z') of the lattice points in the 

bounding boxes BB corresponding to the meshes M10, Mil, 
M15, M16, M20, M21 and M25 are plus (step S32), and 
thus the signed minimum distances D(r— >) are determined 
without re-division of meshes. This is because these 

25 meshes are present on the outside of the object OB to 

be treated. On the other hand, since all the signed 
functions sgn(z f ) of the lattice points in the bounding 



boxes BB corresponding to the meshes M13 and M18 are 
minus (step S33), the signed minimum distances D(r— >) 
are also determined without re-division of meshes. By 
contrast, these meshes are present on the inside of the 
object OB. The bounding boxes BB corresponding to the 
meshes M9, M12, M14, M17, M19, M22, M23 and M24 include 
outline parts of the object OB. Accordingly, like the 
meshes M7 and M8, the re-division of meshes is repeated 
to determine the signed minimum distances D(r— 

After the signed minimum distances D(r~ *) for all 
the bounding boxes BB are found (step S38), the process 
of step S21, which has been described in connection 
with the first embodiment, is executed, and the to-be- 
treated object OB is rendered by the implicit-function 
representation . 

The following advantageous effect (3) can be 
obtained by the image processing method, image proc- 
essing program product and image processing computer 
according to the present embodiment, in addition to 
the advantageous effects (1) and (2) that have been 
described in connection with the first embodiment: 

(3) The object can be represented more precisely, 
while the amount of computations can be reduced. 

According to this embodiment, the size of the 
bounding box is varied on a location-by-location basis. 
In other words, whether the bounding box includes an 
object outline or not is determined on the basis of 



the signed function sgn(z'). As regards the region 
including an outline part of the object, the signed 
minimum distances D(r— *) are found from many lattice 
points using a plurality of small bounding boxes. 
On the other hand, as regards the region including 
no outline part of the object, the signed minimum 
distances D(r— ») are found from a small number of 
lattice points using a large bounding box. In order 
to render the to-be-treated object by the implicit- 
function representation, it should suffice if the 
positions on the surface thereof are recognized. 
Thus, even if the number of lattice points associated 
with the region including no outline part is reduced, 
the precision of the object surface is not adversely 
affected . 

By greatly reducing the number of lattice 
points associated with the region including no outline 
part, the amount of computations in the CPU 11 can be 
reduced. At the same time, by increasing the number 
of lattice points associated with the region including 
an outline part, the to-be-treated object can be 
represented with higher precision. 

Next, referring to FIG. 32, a description is 
given of an image processing method, an image 
processing program product and an image processing 
computer according to a third embodiment of the 
invention. The third embodiment relates to such a 
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case that a plurality of triangular-shaped polygons 
are included in one bounding box in the first 
embodiment and the second embodiment. FIG. 32 is a 
flow chart illustrating a part of an image process 
5 method according to this embodiment. 

As is shown in FIG. 32, the process up to step 
S15 is executed by the method that has been described 
in connection with the first and second embodiments. 
Thereby, the signed minimum distance D(r— ») from any 

10 one of lattice points to a triangular-shaped polygon 

is found. Subsequently, it is determined whether a 
plurality of triangular-shaped polygons TP are included 
in the associated bounding box BB (step S40) . FIG. 33 
shows an example in which six triangular-shaped polygons 

15 TP1 to TP6 are included in one bounding box BB . 

If only one triangular-shaped polygon TP is 
included in the bounding box BB, the process goes 
to step S16. That is, the CPU 11 executes the same 
process as in the first embodiment. On the other hand, 

20 if a plurality of triangular-shaped polygons TP1 to 

TP6 are included in the bounding box BB, as shown in 
FIG. 33, the CPU 11 confirms whether the signed minimum 
distance D(r— >) associated with the present lattice 
point is already found by a triangular-shaped polygon 

25 (e.g. TP2-TP6 ) that is other than a triangular-shaped 

polygon (e.g. TP1) that is currently of interest (step 
S41). If the signed minimum distance D(r— >) is not yet 



found by the other triangular-shaped polygons TP2 to 
TP6, the control advances to step S16 and the process 
for the next lattice point is executed. 

If the signed minimum distance D(r— *) is already 
found by the other triangular-shaped polygons TP2 to 
TP6 (the value in this case is referred to as "previous 
value"), the CPU 11 compares a latest signed minimum 
distance D(r— ») found in association with the present 
triangular-shaped polygon TP1 (this value is referred 
to as "latest value") with the previous value (step 
S42) . The latest value is less than the previous 
value, the CPU 11 discards the previous value and 
updates the signed minimum distance D(r— > ) associated 
with the present lattice point with the latest value. 
The CPU 11 then goes to step S16. 

If the latest value is not less than the previous 
value (step S42), the CPU 11 determines whether the 
latest value is equal to the previous value (step S43) . 
If both are not equal, that is, if the latest value is 
greater than the previous value, the CPU 11 advances to 
step S16 without updating the signed minimum distance 
D(r— >) . In this case, the latest value is discarded. 
If the latest value is equal to the previous value in 
step S43, the CPU 11 confirms a displacement of a 
projection point on the X-Y plane (step S44). 

The step S44 is explained by taking, as 
an example, a case where the signed minimum distance 
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D(r— ») of the triangular-shaped polygon TP1 is already 
known when the distance D(r— ►) of the triangular-shaped 
polygon TP2 is to be found in connection with the 
bounding box BB shown in FIG. 33. If the signed 
5 minimum distance D(r— *) associated with the present 

lattice point is equal between the triangular-shaped 
polygons TP1 and TP2, a displacement between the 
projection point on the X-Y plane and the nearest 
lattice point is compared with respect to each of the 

10 triangular-shaped polygons TP1 and TP2 . As shown in 

FIG. 34, as regards the triangular-shaped polygon TP1 
with the already known signed minimum distance D(r— >), 
a projection point vector from the lattice point LP, 
which is currently of interest, onto the X-Y coordi- 

15 nates is set as rgi 1 ^/ and the nearest point vector 

is set as r n ^— ». In this case, a displacement between 

r 0l'~~* anc * r nl — *r i.e. a distance dl therebetween, is 

Mr n l-> - r 0 i'^l 2 ) 1/2 . 

On the other hand, as shown in FIG. 35, as regards 

20 the triangular-shaped polygon TP2 that is currently of 

interest, a projection point vector from the lattice 
point LP onto the X-Y coordinates is set as rQ2 ,— *' 
and the nearest point vector is set as r n 2~^. In this 
case, a displacement between ro2 f ~~* anc * r n2~ i.e. a 

25 distance d2 therebetween, is ( | r n 2~"* ~~ rg2 ! ~ M 2 ) . 

Assume now that d2 < dl . That is, consider that 
the displacement on the X-Y coordinates is less in 
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connection with the triangular-shaped polygon TP2 that 
is currently of interest than in connection with the 
triangular-shaped polygon TP1 with the known D(r— >). 
In this case, the CPU 11 updates the signed minimum 
5 distance D(r— ») associated with the lattice point 

with the latest value found in connection with the 
triangular-shaped polygon TP2 (step S45) . Then, the 
CPU 11 advances to the process of step S16. On the 
other hand, if d2 > dl, the latest value is discarded 

10 and the signed minimum distance D(r— 0 is not updated. 

The CPU 11 then goes to the process of step S16. 

The following advantageous effect (4) can be 
obtained by the present embodiment, in addition to 
the advantageous effects (1) to (3) that have been 

15 described in connection with the first and second 

embodiments : 

(4) The object can be represented with still 
higher precision. 

In the method according to the present embodiment, 

20 in a case where a position of an object surface is 

given by a plurality of polygons, the surface is repre- 
sented by the implicit function using the less signed 
minimum distance D(r— >). If the signed minimum 
distance D(r— *) is equal between the plural polygons, 

25 the signed minimum distance D(r— ») with a less 

displacement from the projection point on the X-Y plane 
to the triangular-shaped polygon is used. Therefore, 
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the object surface can be represented with higher 
accuracy . 

The present embodiment is particularly effective 
in connection with the second embodiment. In the first 
embodiment, a bounding box is generated for each of 
polygons. In the second embodiment, a bounding box is 
generated for each of predetermined regions, and so it 
is highly possible that a plurality of polygons are 
included in one bounding box. 

Next, a description is given of an image proc- 
essing method and an image processing program product 
according to a fourth embodiment of the invention. In 
this fourth embodiment, the image process described in 
the first to third embodiments is realized by an image 
processing LSI, and not by a software process. FIG. 36 
is a block diagram of a system LSI according to this 
embodiment . 

As is shown in FIG. 36, an image drawing processor 
system LSI 118 comprises a host processor 20, an I/O 
processor 30, a main memory 40 and a graphic processor 
50. The host processor 20 and graphic processor 50 are 
connected to be mutually communicable over a processor 
bus BUS. 

The host processor 20 includes a main processor 
21, I/O units 22, 23 and 24, and a plurality of digital 
signal processors (DSPs) 25. These circuit blocks are 
connected to be mutually communicable over a local 
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network LN1. The main processor 21 controls the opera- 
tions of the respective circuit blocks in the host 
processor 20. The I/O unit 22 executes transmission/ 
reception of data with the outside of the host 
5 processor 20 via the I/O processor 30. The I/O unit 

23 executes data transmission/reception with the main 
memory 40. The I/O unit 24 executes data transmission/ 
reception with the graphic processor 50 via the proc- 
essor bus BUS. The digital signal processors 25 

10 executes signal processing on the basis of data that 

is read in from the main memory 40 or from outside. 

The I/O processor 30 connects the host processor 
20 to, for instance, a general-purpose bus, a periph- 
eral unit such as an HDD or a DVD (Digital Versatile 

15 Disc) drive, and a network. In this case, the HDD 

or DVD drive may be mounted on the LSI 18 or may be 
provided outside the LSI 18. 

The main memory 40 stores a program that is 
necessary for the operation of the host processor 20. 

20 This program is read out of, e.g. an HDD (not shown), 

and loaded in the main memory 40. 

The graphic processor 50 includes a controller 51, 
I/O units 52 and 53, and an arithmetic process unit 54. 
The controller 51 executes, for example, communication 

25 with the host processor 20 and a control of the arith- 

metic process unit 54. The I/O unit 52 controls 
input/output from/to the host processor 20 via the 



42 



processor bus BUS. The I/O unit 53 controls input/ 
output from/to various general-purpose buses such as a 
PCI (Peripheral Component Interconnect), input/output 
of video and audio, and input/output from/to an 
5 external memory. The arithmetic process unit 54 

executes an image processing arithmetic operation. 

The arithmetic process unit 54 includes a raster- 
izer 55, and a plurality of digital signal processors 
(DSPs) 56-0 to 56-31. In this embodiment, the number 

10 of digital signal processors is set at 32, but this 

number is not limited. Alternatively, the number of 
digital signal processors may be 8, 16, 64, etc. The 
detailed structure of the arithmetic process unit 54 is 
described referring to FIG . 37. FIG. 37 is a block 

15 diagram of the graphic .processor 50. 

As is shown in FIG. 37, the arithmetic process 
unit 54 includes the rasterizer 55, and 32 digital 
signal processors 56-0 to 56-31. The rasterizer 55 
generates pixels according to input graphic informa- 

20 tion. The pixel, in this context, is a minimum-unit 

area that is manipulated when a predetermined graphical 
figure is drawn. A figure is drawn by a set of pixels. 
The pixels to be generated are determined by the shape 
of a figure (i.e. positions occupied by the figure). 

25 To be more specific, when drawing is effected at a 

given position, a pixel corresponding to this position 
is generated. When drawing is effected at another 



position, another pixel corresponding to this position 
is generated. The digital signal processors 56-0 to 
56-31 include pixel process units PPUO to PPU31 and 
local memories LMO to LM31, which are associated with 
the pixel process units PPUO to PPU31 , respectively. 

Each of the pixel process units PPUO to PPU31 
includes four realize pipes RP that constitute a single 
RP cluster RPC (realize pipe cluster) . The RP cluster 
RPC executes an SIMD (Single Instruction Multiple Data) 
operation and processes four pixels at a time. Pixels, 
which correspond to respective positions of a graphical 
figure, are assigned to the pixel process units PPUO to 
PPU31. The pixel process units PPUO to PPU31 process 
the associated pixels in accordance with the positions 
of the figure . 

The local memories LMO to LM31 store pixel data 
that are generated by the pixel process units PPUO to 
PPU31. The local memories LMO to LM31, as a whole, 
constitute a realize memory. The realize memory is, 
for instance, a DRAM. The DRAM includes memory areas 
each having a predetermined data width, and these 
memory areas correspond to the local memories LMO 
to LM31. 

In the graphic processor 50 with the above- 
described structure, polygonized object data is input 
to the I/O unit 52 or 53. Based on the input data, the 
rasterizer 55 executes a division process (step Sll, 



512) for dividing an n-angular <n > 3) polygon into 
triangular-shaped polygons, an extraction process (step 

513) for extracting vertex coordinates of a triangular- 
shaped polygon, and a setting process (step S14) for 
setting a bounding box or a bounding sphere. The 
rasterizer 55 delivers the obtained information to the 
digital signal processors 56-0 to 56-31. The digital 
signal processors 56-0 to 56-31 execute parallel proc- 
essors for different bounding boxes (polygons) (steps 
S15 to S18). In the digital signal processors 56-0 to 
56-31, the plural realize pipes RP can execute parallel 
processes for a plurality of lattice points. On the 
basis of the signed minimum distances D(r— *) obtained 
by the digital signal processors 56-0 to 56-31, an 
object is rendered by implicit-function representation. 
The signed minimum distances D(r— >) are stored in the 
local memories LMO to LM31. In addition, the digital 
signal processors 56-0 to 56-31 can execute the process 
of steps S32 to S34, which has been described in 
connection with the second embodiment. In this case, 
for example, the rasterizer 55 executes the process 

of steps S30, S31, S35, S36 and S37, and the lattice- 
setting number-of-times k is stored in any one of the 
memories. Furthermore, the digital signal processors 
56-0 to 56-31 can execute the process of steps S40 to 
S45, which has been described in connection with the 
third embodiment. 



As has been described above, the processes that 
have been described in connection with the first to 
third embodiments can be carried out using the image 
processor according to the present embodiment. 

According to the image process method, image 
process program product and image processing computer 
of the first to fourth embodiments of the invention, 
the bounding box is generated for each of the 
triangular-shaped polygons, and the signed minimum 
distances D(r— >) from the lattice point to the 
triangular-shaped polygon is found. By overlapping 
the signed minimum distances D(r— *) associated with all 
triangular-shaped polygons, the object is represented. 
Specifically, a set of points with D(r— ») = 0 forms a 
surface of an object, and the set of the points meets 
an implicit function <t> (r— >) = 0. Thus, in the method 
according to the embodiments of the invention, the 
implicit function <t> (r— >) = 0 is not directly found. 
However, by searching for the set of D(r~ *) = 0, the 
implicit function <f> (r— ») = 0 is substantially obtained. 
According to the present method, even where an object 
has a complex shape, an implicit-function 
representation of the object can surely be obtained, 
and an implicit-function representation method with 
high versatility in use can be provided. 

In the method according to the embodiments, the 
object surface is defined by discrete lattice points, 
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as shown in FIG. 7. Thus, when the process illustrated 
in FIG. 10 to FIG. 13 is executed, lattice points are 
not necessarily present on the triangular-shaped 
polygon. In such a case, an approximation process is 
5 performed, an a nearest lattice point that is not on 

the triangular-shaped polygon may be used as a point 
that meets D(r— *) = 0. 

Not only the implicit function 0(r— >), but also 
V <f> (r— *) may be found at the same time. Thereby, 
10 even at a position where no lattice point is present, 

4> (r— > ) can be supplemented by a cubic polynomial and 
its value can be found. In this case, V<Mr— *) is 



given by 



15 V+(r) = 



r dif d§ r ~°* ^' ^' A 



K dx dy dz j 



5 2 4> d 2 § d 2 $ d 3 § 
dxdy ' dydz ' dxdz ' dxdydz 



In the above-described embodiments, the implicit 
function <f> (r— >) = 0 is obtained on the basis of the 
signed minimum distances D(r— >) associated with all 
lattice points included in the bounding box. However, 

20 in order to represent an object surface, it should 

suffice if the implicit function $ (r— >) is obtained 
on the basis of the signed minimum distances D(r— >) 
associated with lattice points that are present within 
a distance of less than f from the triangular-shaped 

25 polygon. As regards lattice points at a distance of 

more than e from the triangular-shaped polygon, it 
should suffice if the signed minimum distance D(r— ») 



47 



has an equal signed function sgn(z') and an absolute 
value greater than s . Thus, <}> (r— >) can be found by 
simple extrapolation. In addition, when <t> (r— >) is 
obtained using the lattice points at a distance of 
5 more than e from the triangular-shaped polygon, re- 

initialization may be executed. The re-initialization, 
in this context, refers to correction that is executed 
when <$> (r— ») of the implicit function <S> (r— ») = 0 fails 
to meet the property of the signed minimum distance 
10 function. The re-initialization is executed by 

performing a steady-state calculation for x in the 
following equation : 

§ = s«i> x= o)(i - N) 

ox 

15 where x is a time step. The re-initialization is 

executed to prevent the shape, which is represented 
by the implicit function, from failing to meet the 
property of the signed distance function when the 
shape varies with time. Accordingly, a variable 

20 relating to time is included in the equation. In 

addition, S ( T _ q) is an approximate function of 
the signed function sgn and is given by 

S (()>) = : t « sgn (<|>) 

25 where 8 is a value that is close to 2 e , and this 

equation is established when 5 is sufficiently small. 

In the above-described embodiments, all polygons 
are triangular-shaped polygons. However, needless to 
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say, a rectangular-shaped polygon may be used as it 
is. In this case, the signed minimum distances D(r— *) 
from two triangular polygons that are included in the 
rectangular-shaped polygon may be overlapped. 

The image process program products, image proc- 
essing computers or system LSIs according to the 
first to fourth embodiments are applicable to, e.g. 
game machines, home servers, TVs, mobile information 
terminals, etc. FIG. 38 is a block diagram of a 
digital board included in a digital TV that executes 
the image processing methods according to the first to 
fourth embodiments. The digital board is employed to 
control communication information such as video/audio. 
As is shown in FIG. 38, the digital board 100 comprises 
a front-end unit 110, an image drawing processor system 
120, a digital input unit 130, A/D converters 140 and 
180, a ghost reduction unit 150, a 3D YC separation 
unit 160, a color decoder 170, a LAN process LSI 190, a 
LAN terminal 200, a bridge media controller LSI 210, a 
card slot 220, a flash memory 230, and a large-capacity 
memory (e.g. dynamic random access memory (DRAM)) 
240. The front-end unit 110 includes digital tuner 

« 

modules 111 and 112, an OFDM (Orthogonal Frequency 
Division Multiplex) demodulation unit 113, and a QPSK 
(Quadrature Phase Shift Keying) demodulation unit 114. 

The image drawing processor system 120 comprises a 
transmission/reception circuit 121, an MPEG2 decoder 



122, a graphic engine 123, a digital format converter 
124, and a processor 125. For example, the graphic 
engine 123 and processor 125 correspond to the graphic 
processor 50 and host processor 20, which have been 
described in connection with the first to fourth 
embodiments . 

In the above structure, terrestrial digital 
broadcasting waves, BS (Broadcast Satellite) digital 
broadcasting waves and 110-degree CS (Communications 
Satellite) digital broadcasting waves are demodulated 
by the front-end unit 110. In addition, terrestrial 
analog broadcasting waves and DVD/VTR signals are 
decoded by the 3D YC separation unit 160 and color 
decoder 170. The demodulated/decoded signals are 
input to the image drawing processor system 120 and 
are separated into video, audio and data by the 
transmission/reception circuit 121. As regards the 
video, video information is input to the graphic 
engine 123 via the MPEG2 decoder 122. The graphic 
engine 123 then renders an object by implicit-function 
representation by the method as described in the 
embodiments . 

FIG. 39 is a block diagram of a recording/ 
reproducing apparatus that executes the image 
processing methods as described in connection with the 
first to fourth embodiments. As is shown in FIG. 39, a 
recording/reproducing apparatus 300 comprises a head 
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amplifier 310, a motor driver 320, a memory 330, an 
image information control circuit 340, a user I/F CPU 
350, a flash memory 360, a display 370, a video output 
unit 380, and an audio output unit 390, 

The image information control circuit 340 include 
a memory interface 341, a digital signal processor 342 
a processor 343, an audio processor 344, and a video 
processor 345. 

With the above structure, video data that is 
read out of the head amplifier 310 is input to the 
image information control circuit 340. Then, graphic 
information is input from the digital signal processor 
342 to the video processor 345. The video processor 
345 draws an object by the methods described in the 
embodiments of the invention. 

Additional advantages and modifications will 
readily occur to those skilled in the art. Therefore, 
the invention in its broader aspects is not limited 
to the specific details and representative embodiments 
shown and described herein. Accordingly, various 
modifications may be made without departing from the 
spirit or scope of the general inventive concept as 
defined by the appended claims and their equivalents. 



